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ABSTRACT 


The energization process of magnetic pumping, a combination 
of time-dependent magnetic mirror fields with pitch-angle scatter- 
ing, is applied to trapped charged particles 2x3 drifting in 
corotating, azimuthally nonsymmetric neutron star magnetospheres. 
When particle energization is balanced by synchrotron radiation 
loss, it is found that protons, rather than electrons, reach con- 
siderable kinetic energies and radiate, In the x-ray regime, at 
rates up to the 10^ MeV/proton/sec. 

Sub.iect headings : plasma processes — x-ray binaries — 

magnetospheres — neutron stars 
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I. INTRODUCTION 

Several compact x-ray sources are known to be members of 
close binary systems (Gurskey and Schreir 1975)- We propose that 
the environment of a magnetized neutron star in a close binary sys- 
tem is favorable for the energization, by the process of magnetic 
pumping, of the magnetospheric charged particles, and that these 
particles will be capable of emitting, with considerable power, 
synchrotron radiation in the x-ray region of the spectrum. 

The magnetic-pumping mechanism for the energization of 
charged particles in magnetic-mirror geometries was first proposed 
by Alfv^n (1954) . In his simple model, described in Alfv^n and 
FAlthammar (1963, Sect. 2.7-4), cyclic variations of the magnetic 
field strength are combined with isotropization of the particle dis- 
tribution functions when the field strength reaches a maximum or a 
minimum. As the field strength is increased, the mean perpendicular 
(to 5) momentum of the distribution increases, conserving the first 
adiabatic invariant of each particle, while the mean parallel momen- 
tum remains constant. Isotropization when the field is at a maximum 
turns some of the high perpendicular momenta of the distribution 
into high parallel momenta, which do not decrease with the decreas- 
ing of the magnetic field. This leads to an increase in the mean 
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momentum of the distribution with time. In the Alfvenic model, the 
momentum after one cycle Is given by 



5 + 2k + 2/k 



( 1 ) 


where P Q Is the initial momentum and k is the magnetic ratio of the 
cycle, 


k = B /B 

°MAX °MIN 


(2) 


the B q values being the field strength in the magnetic mirror mid- 
plane (equatorial region for magnetospheres). Equation (l) results 
in a variation in time for the mean momentum, if the magnetic 
pumping is not balanced by another mechanism, which is of the form 




f 


( 3 ) 


where t e is defined as an energization time. 

More realistic, but not simple to solve, is the case of 
pitch-angle (momentum-conserving) scattering of particles occurring 
constantly as the magnetic field varies, in contrast to the infinite 
pitch-angle diffusion imposed at times of magnetic field maxima and 
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minima in the Alfv^n model. In order to solve this problem, we have 
developed formalisms to follow the adiabatic motions of particles in 
the magnetic-mirror topology of an asymmetric magnetosphere and to 
solve a diffusion equation representing the pitch-angle scattering 
of particles. 

As the magnetic field varies we numerically follow, in 
momentum-pitch- angle phase space, the evolution of n(p,a Q ), the num- 
ber of particles in a flux tube with momentum p and equatorial (at 
the mirror midplane) pitch- angle a Q , calculating the new phase- space 
coordinates. In our relativistic treatment, particle motion con- 
serves the first adiabatic invariant 



( 4 ) 


(no longer the magnetic moment), the total energy 


2 

E = 7 me ± eV , 


and the second adiabatic invariant 


( 5 ) 


J = # P„ ds 


t 


( 6 ) 
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where i and \\ are with respect to the magnetic field, V is an elec- 
trie potential, and y is the relativistic factor (l - V*“/c ) ' . 

In order to simplify the treatment, a special shape for the 
magnetic-mirror field is chosen (Goertz 1978; Borovsky et al. 1980), 
with the magnetic -pumping results found to be quite insensitive to 
the choice of magnetic field topology. For a more detailed discus- 
sion we refer the reader to the paper by Borovsky et al. (1980) 
where a similar mechanism was applied to the case of Jupiter. 

In order to describe the diffusion, caused by particle 
scattering in the distribution, in pitch-angle space, we solve the 
"bounce- averaged" Fokker-Flanck equation 

It - £r « 0 “o £r 

o o 

x [n(a o )/T(a o )sin a Q cos a Q )] , (7) 

where T(a Q ) is the angular dependence of a particle bounce period 
and Jfi( a 0 ) is the "bounce- averaged" diffusion coefficient. Bounce 
averaging allows the description of off-equatorial (away from the 
midplane) scattering in terms of equatorial parameters. The time 
asymptotic solution to the diffusion equation is an isotropic dis- 
tribution function f(a Q ) = constant, where f (a^) and n(a Q ) are 


related by 
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f(a ) b n(a )/T(a )sln a cos a 
o o' o o o 


( 8 ) 


Combining the adiabatic motion with the pitch- angle scatter- 
ing allows us to simulate magnetic -pumping processes. For the 
numerical procedure chosen, it is easily shown that alternately 
solving for the adiabatic motion and for the pitch-angle diffusion 
is equivalent to solving the Boltzmann equation for particles under- 
going the magnetic -pumping process. In conforming our simulation 
parameters to the special case of the Alfv4n model, we find very 
satisfactory agreement between our results and analytical theory. 


II. THE MAGNETIC FIELD 


The model magnetic field chosen is a rotationally aligned 
(magnetic moment parallel to angular momentum) dipole with a slight 
compression (increase) on one side. The equatorial field strength 
on the dipolar (noncompressed) side is given by 

■ ( 9 ) 


i 


where M is the magnetic dipole moment. For the form of the 
compressed-dipole side, we choose 



( 10 ) 


A measure of the magnetic field distortion is r CQM p, the radius 
where the compression amounts to a ICffo increase in the equatorial 
field strength, i.e.. 


B = 1.1 at r 
c v 


r C0MP 


( 11 ) 
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yielding, from (10), 


M' 


r C0MP 



( 12 ) 


Thus, the values of M and will define a value for M\ 

We assume that the potentials in the Inner magnetosphere 
arise from corotational electric fields (we discuss corotation 
below) and we will require that the potentials of both sides be 
edaal at the atelier surface, r = V On the dipolar side we have 
(hy Integrating E = l/c v x 1L) , 


y fi M 
''D “ c r 


(15) 


where Si Is the angular frequency of stellar rotation, and on the 
compressed side we have 


V 

c 


fi [m 
c r 


+ M' log 



(14) 


which is equal to Vp at r = r g . 

The adiabatic motion of charged particles conserves (Equa- 
tions (4), (5), and (6)) the first adiabatic invariant p., the 
second adiabatic invariant J, and the total energy E. For particles 
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with much lees kinetic energy than potential energy, the drift 

orbits will lie almost exactly on equlpotentlal surfaces and the 

particles will drift through regions of differing magnetic field 

strength. (The violation of this will provide a limit for the 

energizatlou process, discussed below. ) Since we know the magnetic 

fields B(r) and potentials V(r) on both sides of the star, we may 

calculate the ratio k » B /B for particles in orbits follow- 

°MAX °Min 

ing equlpotentlals. Hence, for choices of diffusion coefficients 

(see below), we may calculate the magnetic-pumping energization 

times, f E , as functions of radius. 

We estimate, for the parameters of our model, the regions of 

the magnetosphere in which the model will be plausible. We take the 

6 

(neutron) stellar radius to be 10 kilometers (r u ~ 10 cm) . This 

o 

will be the inner limit for the validity of the model. We do not 

expect corotation of the magnetosphere out to the light cylinder 

(radius at which the corotational velocity equals the velocity of 

light), which for our rotational period = 2 sec, is 

r L = 95, 500 km = 9*55 X 10 ^ cm. Corotation might be expected (Lamb 

et al. 1975) to hold out to the distance where the corotational 

velocity is equal to the Kepplarian orbital velocity (r > r^ "» 

> V-^-rm) • This Kepplarlon radius is given by 
CUR OKB11 



XI 


where m s lo the stellar mass. For m s * 1 we find » 3480 km 
and for m B « .1 m Q we find r^ ■ 1613 km. A relativist ic plasma 
theory for pulsar magnetospheres (Hints and Jackt^n 1974) predicts 
the outer limit of corotation to he 


r COR 



m T 


3/5 r 2/5 
s r L 


f 


(16) 
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which, for our model, is r rQR * 391 km 3 3> 91 X 10 cm. At these 
limiting radii the corotationiil velocities are nonrelatlvistic 
(r = 391 km * ^roR “ an^t r 3 300 km ■» Vqq R ■ e/320 ), so we may 

assume the particle drifts to he adiabatic, and in particular, may 
be confident that the second adiabatic invariant, J « j> P ds, is 
conserved. 

As an agent for the pitch-angle scattering of the charged 
particles we assume the presence of plasma waves, excited by 
anisotropic particle distributions produced by the adiabatic motion. 
In calculating the energization times as functions of the equatorial 
distance from the star (k(r)known), we use the continuous diffusion 
model and assume that the mechanism runs at its optimum diffusion 
rate ( t dif / t rot = where t dif is ' bhe diffusion time and t Rqt is 
the rotation period, see Borovsky, Chap. V; Borovsky et al. 1980) . 
This choice is not critical since diffusion rates an order of magni- 
tude above or below this value give energization times differing by 
at most, a factor of three. 
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III. RESULTS OF THE MODEL 

We now have the Information necessary for the construction of 
a model which energizes particles in the asymmetric magnetosphere. 
Without additional processes we would expect the energy of particles 
to Increase indefinitely. There are, however# two phenomena which 
limit particle energy. The first is the condition that a particle 
have less kinetic energy than potential energy. If this requirement 
is met# a particle's drift orbit will approximate an electrical 
equipotentl al and the drift will carry the particle through regions 
of differing magnetic field strength# essential for the success of 
the energization process. If# Instead# a particle's kinetic energy 
exceeds its potential energy# then the drift path will follow con- 
tours of equal magnetic field strength and no energization will 
occur (k -♦ l) . To estimate this kinetic-energy upper-limit, we 
calculate the difference in electrical potential# AV# along a con- 
tour of equal equatorial magnetic-field strength and require parti- 
cle kinetic energies to be less than e AV. Thus our limit appears 
as 


" l)mc = e AV 


> 


( 17 ) 
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where 7 is the relativistic factor (l - v S /c 2 )’ 1 / 2 and the l sub- 
script pertains to the value at the limit. This limiting kinetic 
energy appears (in terms of 7j) in Figures 2 - 6 as the hollow 
circles . 

But even when 7 < 7^ the energy does not necessarily Increase 
because energy is lost through synchrotron radiation, the power 
radiated being approximately proportional to the square of the par- 
ticle energy (P « |> 2 - 1] ) . The fact that, in the synchrotron 
mode, the power radiated by a particle of set energy is proportional 
to the inverse- four power of the particle mass, and thus more effec- 
tive for electrons than protons, leads to an Important prediction of 
our model; that protons, not electrons, gain the energy and radiate 
it away. To find the limiting kinetic energy we will equate the 
energization time, Tg, with the radiative lifetime for relativistic 
particles (Alfven and Faith ammar 1963) 



, 3 3 
3 a m 



X 1 

B 2 1 + 7 r 


( 18 ) 


Equating this to Tg yields the relativistic factor 7 as a function 
of the radius r (on the compressed side), 


R 


3 c^ rtP i 1_ 


2 e 


B 2 T E 


- 1 


> 


(19) 
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where B and t e are functions of r. ThlB value appears as the solid 

circles in Figures 2-6. If, at a particular radius, the value 

obtained by balancing energization with radiative loss is lower than 

the value of the mechanism limit (kinetic less than potential), 

then 7 ^ will be the relativistic factor of the particle there, that 

is, particles are energized up to and maintained at the energy 
2 

E a me . If the value is lower than the value 7 R , the partl- 

2 

cles are energized to and maintained at the energy E = y T me . The 
variation of y is shown as the solid line in Figures 2-6, where y 
is always the smaller of the two values, and y^. 

Knowing the particle kinetic energies, (y - l)mc , and mag- 
netic field strength, B, as functions of r, we may calculate the 
wavelength, \ m , of the maximum synchrotron-radiated power as a func- 
tion of r. With the frequency of the peak power being 



el 2 
me ! 


( 20 ) 


we have (\ = c/v) 


2t r mc^ 1 _1_ 
e B 2 

y 


(21) 


These wavelengths, for protons, appear in Figures 2 - 6 as solid 
squares for the radiation-limited proton energies (^ m (?p)) and as 
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hollow squares for the mechanism-limited proton energies (\ m ( 7 L ))> 
the solid curve taking the wavelengths pertaining to the charac- 
teristic (lesser of 7 ^ and 7 ^) value. Electron radiation-limited 
wavelengths appear in these graphs as the solid triangles, hut, due 
to the fact that rapid radiative loss limits them to low kinetic 
energies, the emission is strictly through the cyclotron mode. 

We examine five cases (Figures 2-6), combinations of medium 

28 3 27 3 

(M = 5 X 10 gauss cnr ) or weak (M = 5 X 10 gauss cnr) magnetic 

dipole moments with weak ( r C0Mp = 955 km) or strong (r COMp = 

2 q 3 

95-5 km) compressions and one of a strong (M = 5 X 10 ~ gauss cm ) 
magnetic dipole with weak compression. In the weak dipole (meaning, 
also, weak corotational potentials) cases (Figures 2 and 5), the 
proton energies are more restricted by the (kinetic less than poten- 
tial) mechanism limit, while, in the medium dipole cases (Figures 3 
and 6 ), protons are energized at distances further from the star, in 
regions in which we are less confident in the application of our 
model. In the case of the medium dipole-weak compression (Fig- 
ure 3 ), the energy limit dominates the radiative balance at r = 

320 km (still within the Hinata and Jackson corotation limit, r CQR ), 
where the characteristic wavelength of the synchrotron emission is 
2.8 X 10 ' cm, in the weak x-ray region. For the case of the medium 
dipole-strong compression (Figure 6 ), mechanism-limiting occurs at 
a distance of r = 234 km, where the characteristic radiation is in 

-Q 

the strong x-ray regime, \ = 1.1 X 10 7 cm. In a strong 
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29 3 

dipole- (M = 5 X 10 gauss cm ; applicable to the Crab Nebula 
pulsar) weak compression (again, r C0Mp . 955 km) Casa (Figure k), we 
find a limited wavelength of 2.2 X 10 ^ cm (strong x-ray region) at 
a distance r = 730 km, which, however, may not be in the corotating 
zone. 

As a value for the power radiated per particle, we use the 
particle kinetic energy divided by its radiative lifetime, 

p _ (l “ A l rcp .,, f (22) 

t R 

where 7 is the appropriate relativistic factor, 7 ^ or 7^, and the 
lifetime is given by Equation (18) (using 7 R or 7^) . The power 
radiated per proton is displayed, for the four above cases, in 
Figure 7 , where the substantial radiated-power prediction 
(> lcA MeV/sec/proton) of the model may be seen. It is important to 
realize that, in the cases studied, the power radiated by electrons 
is less by 10 or more orders of magnitude. This is due to the short 
radiative lifetimes of electrons which keep their kinetic energies 
(and momenta) low. The energy gain per unit time, and hence the 
energy radiated per unit time, is proportional to the particle 
kinetic energy (see Equation (3)), so the low-energy electrons gain 
little energy from the magnetic pumping, hence radiate little power. 
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We note that (for 7 » l) as the stellar rotation rate 
ft = 2 tt T rqt Increases, the radiation-limi'Ufl 7- values, y R (Equa- 
tion (19) ) , increase proportional to ft, the characteristic synchro- 
tron wavelengths, (Equation (21)), decrease proportional to ft’ 2 , 

and the power radiated per particle, P (Equation (22)), Increases 
2 

proportional to n . We also note that, since the corotational elec- 
trical potentials (Equation (15)) are proportional to the angular 
rotation rate ft, the mechanism- limiting 7- values, 7 L (Equation (17 ) ) , 
increase proportional to ft as do the radiation- limited values, y R . 
Thus, if we increase the rotation rate of our model by a factor of 
n/lO (ft ■* 1JI0 ft and T R01 , = 2 sec -* t rqt = .65 sec), the characteris- 
tic wavelengths as functions of the radius will all decrease by a 
factor of 10 and the radiated power per particle will increase by a 


factor of 10. 
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IV. DISCUSSION 

A numerical model of the magnetic -pumping process was con- 
structed "by alternately following the adiabatic motion of particles 
trapped in magnetic mirrors and solving a pitch-angle diffusion 
equation. The results of simulations suitably agreed with theo- 
retical predictions for the special case of the Alfv4nic model. 

Our numerical techniques allowed us to manage the magnetic- 
pumping process when pitch-angle diffusion occurs simultaneously 
with the magnetic field variation. For this case we find that, as 
in the Alfvenic model, the mean momentum of a distribution of 
trapped particles, on the average, increased with time t according 
to 

P(t) = P(t = 0)e t//T , (23) 

where r is a positive constant. A resonance between the magnetic- 
field cycle and the pitch- angle diffusion rate gives a maximum 
energization efficiency for the mechanism, that is, for a certain 
value of the (diffusion per magnetic cycle (regulated by a diffusion 
coefficient), the constant t will be a minimum. 
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Rotating neutron stars with slightly asymmetric magneto- 
spheres were modeled. By balancing the energy gain from the 
magnetic -pumping process with the energy loss by synchrotron radia- 
tion, we determined the energies, synchrotron wavelengths and power 
outputs of particles in the inner magnetospheres. We expect syn- 
chrotron radiation by protons, rather than by electrons, and for 

the parameters of our models, predict energy outputs in the range of 

2 6 

10 - 10 MeV per second per proton. 

In order to estimate the luminosity of the neutron star, we 
must estimate the particle density in its magnetosphere. For an 
upper limit to this density we require that the plasma diamagnetic 
field be less than the magneto spheric field (P s 1 in plasma 
terminology), written 

particle pressure s magnetic pressure , (24) 

or, since the electron pressure is negligible, 

2 B 2 

n(r “ l)mc s: ^ , (25) 

where m is the protonic mass and B is the magnetospheric field. 

This gives an upper limit for the particle number density of 
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8n(y - l)mc' 


2 * 


This limit is used to estimate the magnetospheric luminosity 
interior to the distance r, defined 


L(r) = J P(r)n(r) dr 


where the power radiated per particle, P(r), is obtained from Equa- 
tion (22), the factor of l/2 comes from assuming half the spherical 
volume contains particles, and the lower limit r g is the neutron- 
stellar radius. We display this luminosity L(r) for various field 
strengths, etc., in Figure 8. Note that since the plasma diamag- 
netic effect density limit (26) is independent of the angular rota- 
tion ratio fi, and since this limit for n is several orders of 
magnitude smaller than the force balance limit given by (in = mass 
neutron star) 


mm /R 2 n 2 ? 

G ~ 2 ~ n - ^q"^/ “ ?(n(y - l)mc ) = nmfl r , 
r 


j 

/ 
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an Increase in fl by a factor of <JlO will change P(r) by a factor of 
10 without changing n(r), hence will Increase the luminosity L(r) by 
a factor of 10. 

Finally we would like to point out that other losses besides 
synchrotron radiation (e.g., outward mass transport, collisions! 
losses, precipitation onto the star) have not been considered in our 
model. All these losses are only important if the loss times asso- 
ciated with the processes are shorter than the energization time or 
radiation lifetime. Collision times are certainly much larger than 
the times considered here which are of the order of several rotation 
periods. Only strong mass transport rates, such that the plasma 
radial velocity is larger than the corotation speed, would reduce 
the efficiency of our mechanism. The effect of precipitation is 
difficult to estimate but the extreme smallness of the loss cone in 
the region where the synchrotron radiation is generated should 
guarantee a small loss rate. 
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FIGURE CAPTIONS 

Figure 1. Time dependence of mean momentum of a distribution In the 
continuous-diffusion magnetic-pumping model. 

Figure 2. Neutron star magnetospherlc model: weak magnetic dipole 

moment, weak compression. 

Figure 3. Neutron star magnetospherlc model: medium magnetic 

dipole moment, weak compression. 

Figure 4. Neutron star magnetospherlc model: strong magnetic 

dipole moment (applicable to Crab Pulsar), weak compression. 

Figure 5* Neutron star magnetospherlc model: weak magnetic dipole 

moment, strong compression. 

Figure 6 . Neutron star magnetospherlc model: medium magnetic 

dipole moment, strong compression. 

Figure 7* Power radiated per proton in the synchrotron mode for 
various magnetospherlc parameters. 

Figure 8. Synchrotron luminosities of magnetospheres interior to 
radius r as functions of r. 
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